Theoretical Study of Nuclear Spin Polarization and Depolarization in Self- Assembled 

Quantum Dots 
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We investigate how the strain-induced nuclear quadrupole interaction influences the degree of 
nuclear spin polarization in self-assembled quantum dots. Our calculation shows that the achievable 
nuclear spin polarization in In^Gai-^As quantum dots is related to the concentration of indium 
and the resulting strain distribution in the dots. The interplay between the nuclear quadrupole 
interaction and Zeeman splitting leads to interesting features in the magnetic field dependence of 
the nuclear spin polarization. Our results are in qualitative agreement with measured nuclear spin 
polarization by various experimental groups. 



I. INTRODUCTION 

Nuclear spin dynamics has been studied extensively 
in many fields In recent years, nuclear spin dynamics 
in semiconductor quantum dots has attracted intense in- 
terest because of the excellent quantum coherence prop- 
erties of nuclear spins. Indeed, nuclear spins in nanos- 
tructures have been suggested as qubits for a quantum 
computer 3 and for use as quantum memory^ For either 
purpose, high degree of nuclear spin polarization (NSP) 
is a pre-requisite. 

Nuclear spins in nanostructures also form an impor- 
tant environment for electron spins, which have been 
proposed candidate for qubits £ Through the hy- 

perfine interaction the nuclear spins create a spatially 
and temporally fluctuating magnetic field for the elec- 
tron spins, which leads to spin decoherence£~— It has 
been suggested theoretically that such decoherence could 
potentially be suppressed if nuclear spin fluctuations are 
suppressed^ - — and one way to realize such suppression 
is via dynamic nuclear spin polarization. Furthermore, 
the coupled electron-nuclear spin problem is an intrigu- 
ing example of a quantum many-body problem, and is 
still not solved completely. 

Dynamic nuclear spin polarization (DNSP) has been 
studied for many decades. 15 It has been demonstrated 
in semiconductor quantum wellsi^— and quantum 
dots^2r— through a variety of experiments. Physically, 
DNSP can be achieved either electrically or optically, 
where a pumped electron can transfer its spin polariza- 
tion to nuclear spins via the contact hyperfine interaction. 
A range of values for nuclear spin polarization has been 
reported by several experimental groups. For example, 
using electrically controlled DNSP, Petta et al. reported 
approximately 1 % NSP in lateral coupled GaAs dou- 
ble quantum dots^ 3 . while Baugh et al. reported 40 % 
NSP in vertical coupled GaAs quantum dots with 5% 
InjSi With optically pumped DNSP, Gammon et al. re- 
ported 60 % NSP in interface fluctuation GaAs quan- 
tum dotsji 9 - while recent experiments have achieved NSP 
in In x Gai_ x As self-assembled quantum dots at various 
magnetic fields. In particular, approximately 10 ~ 20% 
of NSP is created in Ino.9Gao.1As quantum dots < 1 
T^ 40% in In0.eGa0.4As at around 2 T^rM. and 80% 



in Ino.9Gao.1As at 5 T. 28 It is evident that the exper- 
imental results vary greatly as experimental conditions 
and physical systems are varied. So far there has been 
no systematic theoretical studies of NSP and how it de- 
pends on the various system parameters such as applied 
field and material composition. 

In this paper we study dynamic nuclear spin polar- 
ization in In-rGai-^As quantum dots via optical pump- 
ing of confined electrons. These self-assembled dots are 
generally highly strained, and we are particularly inter- 
ested in the NSP of these dots in different strain environ- 
ments. Specifically, the strain breaks the cubic symmetry 
of the crystal lattice and creates an electric field gradi- 
ent which couples to the nuclear quadrupole moment, 29 
which in turn leads to mixing of nuclear spin eigen- 
states. We use a simplified model of the quantum dot 
where the electric field gradient is axially symmetric. We 
first study how the As NSP responds to various strain 
strengths, angles and cotunneling constants in magnetic 
fields. For the electron-nuclear spin transfer, we con- 
sider both phonon-assisted and cotunneling-assisted spin 
flip processes. Lastly, we consider NSP of In^Gai-^As 
quantum dots with different compositions. 

The paper is organized as follows. We describe the 
scheme of DNSP in section HH and our model Hamilto- 
nian in section Mil We show our results of nuclear spin 
polarization in As nuclei and in different compositions of 
In^Gai-^As quantum dots at various magnetic fields in 
section HVl We discuss some interesting features related 
to our calculation in section |Vl Finally, we summarize 
our results in section IVll and draw our conclusions. 



II. SCHEME OF NUCLEAR SPIN 
POLARIZATION 

In r Gai_ x As self-assembled quantum dots (SAQDs) 
are formed by a strain-driven process, where the strain 
arises from the lattice mismatch between the InAs depo- 
sition layers and the GaAs substrate. The strain in the 
quantum dots (QDs) breaks the lattice symmetry and 
creates electric field gradients in the dots. The shape 
of In x Gai_a;As QDs varies among experiments, ranging 
from pancake-like to pyramid-like and dome-like. The 
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Figure 1. A model of a self-assembled quantum dot. We as- 
sume the field gradients in the self-assembled quantum dot are 
axially symmetric. The largest electric field gradient compo- 
nent Vzz is along the principal axis Z, which deviates from 
the growth direction (the z-axis) by an angle 9. The external 
magnetic field is assumed to be along the z-axis. 
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Figure 2. Schematic of the setup for a typical optical ori- 
entation experiment . 31 ' 32 The self-assembled quantum dot is 
embedded in a Schottky diode heterostructure, where the gate 
voltage is tuned to allow only one electron orbital state below 
the Fermi level. Nuclear spins are polarized by an optically 
pumped electron captured in the quantum dot. 

resulting distribution of electric field gradients thus also 
differs from dot to dot. Even in the same dot, the strain 
distribution is not uniform. For example, the strain at 
the edge of a quantum dot is generally larger than at the 
center of the dot<2£ Therefore nuclear spins in different 
regions of a QD experience electric field gradients of dif- 
ferent strengths and directions. For an estimate of NSP, 
we start with a simplified model of a pancake-like cylin- 
drically symmetric QD as shown in Fig. [T] The electric 
field gradients in such a dot are thus axially symmetric. 
The largest electric field gradient Vzz is along the princi- 
pal axis Z, which is defined to be normal to the pancake 
surface. For instance, for a lattice site in the xz plane, Z 
would be in the xz plane as well, and deviates from the 
growth direction (the z-axis) by an angle of 6. 



We base our calculation on the experimental conditions 
in Refs. M, M, M, HI (see Fig. 0, where the SAQD is 
embedded in a Schottky diode heterostructure, so that 
the charged states of the QD can be controlled. In addi- 
tion, the gate voltage can be tuned to allow zero or one 
charge (electron or hole) on the dot. 

In such a system dynamic nuclear spin polarization is 
realized via optical pumping. A circularly polarized pho- 
ton creates an electron-hole pair, which is then captured 
in the QD as a negative trion ; 22 ' 32 neutral exciton^ or 
positive trioni2£~— When the QD contains one electron, 
this electron is likely spin polarized and can polarize a nu- 
clear spin through the hyperfine interaction. The prob- 
ability to realize this spin transfer process depends on 
the experimental conditions, such as the type of excita- 
tion used (pulsed or CW) and the initial charged states 
in the quantum dot (X°, A"~or X + ). Therefore, an ex- 
perimentally determined factor f e is used to modify the 
electron- nuclear spin transfer probability (see IIIII C). 

In a finite magnetic field the hyperfine-mediated trans- 
fer of polarization from the electron spin to the nuclear 
spins has to be assisted by another process because of 
the large mismatch of the electron and nuclear Zeeman 
energies. In our case we consider phonon-assisted and 
tunneling-assisted processes. More specifically, in the 
Schottky-diode configuration the confined electron can 
flip its spin via cotunneling to the external reservoir (For 
details see IIII C 1|) . This cotunneling-assisted spin flip 
process is efficient at low magnetic fields. In the high 
field regions (or in an isolated QD), the electron-phonon 
interaction provides the more efficient channel to com- 
pensate for the energy mismatch between electron and 
nuclear spins (For details see IIII C 2|) . 

III. MODEL HAMILTONIAN 

The total Hamiltonian for the nuclear spin polarization 
scheme we consider is given as follows, 

Ht = H n + H e + Hhf , 

where H n is the Hamiltonian for the nuclear spins in the 
quantum dot, H e is for the electron spin, and Hhf is the 
hyperfine interaction between the electron and nuclear 
spins. Below we describe each of the terms in Ht in de- 
tail, and discuss the role they play in the DNSP process. 

A. The Hamiltonian of nuclear spin in quantum 
dots 

In the presence of an external magnetic field (along the 
z-axis) , the Hamiltonian of nuclear spins in our simplified 
model of a quantum dot (see Fig.Q} is defined as follows, 

H n = H z n + H Q + H d , (I) 
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Table I. Material parameters used in our calculation. Nuclear 
electric quadrupole moment Q and constant Sn (which re- 
lates electric field gradient to strain) are taken from Ref. |33j 
and Ref.dU. Nuclear spin gyromagnetic ratios are taken from 
Ref. HI 



Elements 


In 


Ga 


As 


nuclear spin / 


9/2 


3/2 


3/2 


electric quadrupole moment Q (10~ 24 cm 2 } 


0.86 


0.27 


0.2 


Sn (10 15 statcoulombs/cm 3 ) 


16.7 


9.1 


13 


gyromagnetic ratio 7 (/neV/T) 


0.039 


0.042 


0.03 


hyperfine constant A (/xeV) 


56 


42 


46 
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H™ represents the nuclear Zeeman energy (E^) where 
I z is the projection of a nuclear spin along the external 
magnetic field, and 7 is the nuclear gyromagnetic ratio 
(see Table Q}. 

Hq represents the electric quadrupole interaction 
(Eq)^^ through which the nuclear spins in an 
In-rGai-^As SAQD couple to the electric field gradi- 
ents in the crystal lattice. The asymmetric part of the 
quadrupole interaction is neglected here because we as- 
sume a pancake-shaped QD, where Vzz 3> Vxx ,Vyy ■ 
Q is the electric quadrupole moment of a nucleus, and e 
is the elementary charge. Vzz is the electric field gradi- 
ent along the principal axis Z. Vzz — S%\ezzr* where 
the constant Sn is experimentally determined (see in Ta- 
ble HI, and ezz is the Z component of the strain tensor, 
which is approximately 4% to 8% for In^Gai^^As QDs 



with different compositions 



30.36-40 



The electric field gra- 



dients introduced by charged states in the QEe 1 - are at 
least two orders of magnitude smaller than the electric 
field gradient caused by the broken symmetry of the crys- 
tal lattice, therefore we exclude the effect of the former. 

Hd represents the dipolar interaction between differ- 
ent nuclear spins, where Rjk is the distance between the 
ith and the jth the nucleus, /zo is the free space perme- 
ability. The local field Bi oc that this dipolar coupling 
produces is around 1 Gauss, and the dynamical effect of 
the dipolar interaction is nuclear spin diffusion. In our 
calculation, we do not deal with the case when the ex- 
ternal magnetic field is smaller than the local field (the 
smallest external field we use is 10 mT). Furthermore, 
nuclear spin diffusion is strongly suppressed in a small 
QDj 4 ^ and its time scale (tens of seconds to minutes 43 ) 
is much longer than the time scale of DNSP in our study. 
Therefore we exclude Hd in our calculation. Since Hd is 
the only direct interaction between nuclear spins, its re- 
moval significantly simplifies our study: the nuclear spins 
can now be treated independently from each other. The 
hyperfine-mediated nuclear spin interactio n 9 ' 44 ' 45 is also 
neglected in our calculation, as the hyperfine interaction 
is only turned on randomly for a small fraction of time 
in experiments, as we will discuss in Section III.C. 

The nuclear spin Hamiltonian is thus simplified as 
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In our simplified model of the quantum dot, the prin- 
cipal axis Z for the largest electric field gradient V Z z 
deviates from the z-axis by an angle 9, and I z is the pro- 
jection of a nuclear spin along the principal axis Z . For a 
pancake-shaped QD, the angle 9 is generally quite small. 
Therefore, while for all the calculations presented in this 
paper we treat the nuclear spin Hamiltonian ([1} exactly, 
for the qualitative discussion in Section IIVI we take a 
small-angle approximation and simplify the Hamiltonian 
((T|). For example, in the case of / = 3/2 we obtain, 
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where E' Q = ,Eq(3cos 2 (9 - l)/2 w E Q . When the off- 
diagonal terms in Eq. ^ are small compared to the 



diagonal terms, we can construct the new eigenstates per- 
turbatively. In a non-degenerate case, the nuclear spin 
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Figure 3. Sketches of the energy-level diagram of a sin- 
gle nuclear spin (not to scale). Due to the presence of 
the quadrupole interaction, the eigenstates are generally no 
longer the eigenvectors of I z . We label the real eigenstates 
as {|1), |2), |3), |4)}, which are combinations of I z eigen- 
states with m z = {§, !,-§) -§}• (a) Low field situations, 
when the quadrupole interaction (Eq) is stronger than the 
nuclear Zeeman effect (E™). (b) 2-4 degeneracy, when the 
quadrupole interaction and nuclear Zeeman energy are res- 
onant (E z = Eq). (c) 3-4 degeneracy, when E™ = 2Eq. 
(d) High magnetic field region, when the Zeeman energy is 
dominant, and the quadrupole interaction is negligible. 



+ Y, Vk ( C L C -+ C Uka) , (4) 
k<7 

H w = X! M i"i a -qu + a qu )exp(i(\ ■ r). 

qv 

Here iT| represents the electron Zeeman energy, where 
g* is the electron effective g-factor in In^Gai-^As QDs, 
fiB is the Bohr magneton, and S z is the z compo- 
nent of the electron spin operator. Hx is the Anderson 
Hamiltonianj 4 ^— suitable for describing the experimen- 
tal setup in our consideration (see Fig. [5]), where a QD 
is tunnel-coupled to the outside Fermi sea. In Eq. 
and c represent electron creation and annihilation opera- 
tors, and n the number operators. We describe a reservoir 
state with index k, energy £/., and electron spin index a. 
The single electron energy level in the quantum dot is e . 
U is the on-site Coulomb interaction and 14 is the tunnel- 
ing matrix element. H ep represents the electron-phonon 
coupling, where a_ qv and a qv represent phonon creation 
and annihilation operators, with quasi-momentum q and 
branch index v. We consider both the deformation po- 
tential and piezoelectric potential in M g „4£ 



eigenstates are as follows: 
3 V 



-d 



(3) 



where a = V3E Q 6/(E^ + 2E Q ), b = V3E Q 6 2 /4(E% + 
E Q ), c = V3E Q 8 2 /4(E? - E Q ) and d = V3E Q 6/(E? - 
2Eq). For a degenerate case as shown in Fig. [3] (b) and 
(c), the eigenstates have to be solved by directly diago- 
nalizing the Hamiltonian in Eq. @ . Taking the state 2-4 
degeneracy, for example, one would expect a complete 
mixing between state 2 and state 4. 



B. The Hamiltonian of electron spin in quantum 
dots 

The effective Hamiltonian of the electron in the quan- 
tum dot consists of three main parts: Zeeman split- 
ting, tunnel coupling to the external reservoir, and the 
electron-phonon interaction: 



H e = H* + Hj 



C. The hyperfine interaction 

In our scheme, nuclear spin polarization (NSP) is 
pumped by optically oriented electrons via the contact 
hyperfine interaction: 

N 



H hf {t) = h{t)Y,M${Rk 



I k z S z + - {I k + S- 



I k S- 



(5) 



where 



where A is the hyperfine coupling constant (see Table fl}. 
N is approximately 10 4 in an In^Gai-^As SAQD. Q(Rk) 
is the electron wave function at the A:th nucleus site, 
which is Gaussian for harmonic confinement. For our 
calculations presented in this paper, we take |$(i?fc)| 2 as 
1/N, effectively assuming a constant electron wave func- 
tion in the QD. This assumption makes the definition of 
nuclear spin polarization well-defined, while in the case of 
a Gaussian wave function the calculation of overall NSP 
depends on where the dot is truncated, as the edge of the 
dot would generally be only slightly polarized, hit) is a 
random function depending on the experimental proce- 
dures and conditions. While we are not going to describe 
the full details of each pumping scheme, the nature of h(t) 
depends on whether the system is in the trion or neutral 
exciton regime. For example, for the X + scheme, h(t) 
is dependent on the trapping and recombination of the 
electron^ while for X~ scheme it is dependent on the 
exciton recombination and electron tunneling time j 32 i 49 
and for X° it is dependent on the exciton recombination 



H z e = -g*^ B B z S z 
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Figure 4. (color online). Comparison of cotunneling-assisted 
nuclear spin pumping rates with phonon-assisted rates. Here 
the Overhauser field is anti-parallel to the external magnetic 
field, and the angle between the largest electric field gradient 
and the external magnetic field is 2° (we use these condi- 
tions in the following figures, unless otherwise noted). The 
cotunneling-assisted spin flip process is more efficient at low 
to intermediate magnetic fields (B < 5 T). For higher fields, 
the phonon-assisted spin flip process is more efficient. 



timej*i The fraction f e is defined as the mean value of 
this temporal function h(t), and it represents the fraction 
of the time when only one electron is left in the quantum 
dot and the hyperfine interaction is "turned on", so that 
the electron- nuclear spin flip-flop can be realized. We 
use 0.035 for f e in our calculations, based on experimen- 
tal observations^ The small value of f e , together with 
the fact that h(t) is random in time to a degree, justify 
our approximation of neglecting higher-order effects of 
the hyperfine interaction throughout our calculations. 

Depending on the helicity of the optical excitation 
((7 + /(7~) relative to the applied external magnetic field, 
and through the flip-flop terms in the hyperfine inter- 
action, the electron can pump nuclear spins either to 
the highest-energy spin state or the lowest-energy state. 
These spin flip-flops are responsible for pumping the nu- 
clear spins in the NSP process (blue one-way arrows in 
Fig.[SJ). In subsection IIII C II and MI C~2l we calculate the 
pumping rates based on the specific physical processes 
involved. 



1. Cotunneling-assisted spin flip processes 

The electron spin in the QD can interact with an elec- 
tron spin in the Fermi sea via cotunneling processes, so 
that the electron Zeeman levels are broadened. The 
spin flip probability and level broadening can be cal- 
culated applying Schrieffer-Wolf transformation to Eq. 
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Figure 5. (color online). Transitions induced by the hyperfine 
interaction among nuclear spin states with / = 3/2. The blue 
one-way arrows represent the transitions induced by the flip- 
flop terms in Eq. (JSJ , and are responsible for pumping of the 
nuclear spins by the trapped electron spin. For an open-dot 
system as shown in Fig. [2] the cotunneling-assisted spin flip 
rate D is dominant at low to intermediate magnetic fields. At 
high fields or in an isolated dot, the phonon-assisted pump- 
ing rate T ph is dominant. The red two-way arrows represent 
transitions due to the mixture of different nuclear spin states, 
which originate from the I Z S Z term in Eq. ((5j|. These pro- 
cesses are responsible for depolarization. 



46.47,50,51 



For each Zeeman level, 



1 



2rr ( £l -e ) 2 + r2/4' 



where r c is the level broadening due to the cotunneling 
processes. 

The probability of electron nuclear spin flip-flop pro- 
cesses increases when the overlap between the two broad- 
ened electron Zeeman states increases. Therefore the 
cotunneling-assisted spin flip-flop is more important at 
low magnetic fields. The transition rate between an ini- 
tial state |t i) and a final state \\. j) can be calculated 
with the Fermi Golden Rule as. 



2tt 



A 



p(ei)p(£f)5(e f - Ei~ A E)d£ide f 
2/ e r c |qj|/ + 5_|t»)| 2 



(6) 



where S n — 2A (I z ) stands for the Overhauser shift, and 
the + sign in front of the Overhauser field is due to 
our choice here that the external magnetic field is anti- 
parallel to the Overhauser field. For parallel fields, — S n 
should be used. The correlation time, r c = l/r c , de- 
scribes the broadening of QD electron states due to co- 
tunnueling processes. It is estimated for a typical Schot- 
tky structure to be ~ 20 ns4£ At low to intermediate 
magnetic fields, the cotunneling-assisted spin flip pro- 
cesses are the most efficient in building up NSP in the 
QD, as shown in Fig. 5] 
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2. Phonon- assisted spin t 



processes 



For an isolated dot, or a dot described in MI C II in 
higher field regions, the cotunneling-assisted spin flip pro- 
cesses become less efficient due to the larger electronic 
Zeeman splitting. Now the phonon-assisted spin flip pro- 
cesses give the most efficient DNSP channel. The pump- 
ing rates due to the phonon-assisted spin flip processes 
are, 



i,0 



2tt x - . 



qv 

x [n qu 5(hsq ~ E%) + {n qu + 1)5(E% + hsq)} 

= N q [f def (El) + f p iezo(E%)] ■ 



where 



T ep — 



{mjj\H hf \ni) (l\H ep \m) 
E m -E l + El 



(m\H ep \l) (nj\H hf \mti) 



E n 



Ei-E% 



m and I Fock-Darwin orbital states. For simplicity, we 
only consider the phonon emission/absorption between 
s orbital and p orbital states (we do not anticipate the 
inclusion of contributions from higher orbital states to 
qualitatively alter our results). 5 is the deformation po- 
tential constant. The piezoelectric constant is denoted 
as eei4 — 2 x 10~ 10 J/m, sound speed: s — 3 x 10 3 
m/s, electron density in InAs: p — 5.7 x 10 3 Kg/m 3 . In 
the presence of an external magnetic field, Fock-Darwin 
energy levels can be represented as Ml± — Ml ± Hu) c /2, 
where Ml is the electronic confinement in the QD and is 
about 30 meV in the type of QD we consider. The cy- 
clotron frequency is lj c = eB/m* where m* = 0.023mo is 
the effective electron mass in the InAs QDs. Iq is the lat- 
eral dimension of the QD. The transition rate due to the 
deformation potential is proportional to the fifth power 
of the electronic Zeeman splitting, while the contribution 
from the piezoelectric interaction is proportional to the 
third power of the electronic Zeeman spitting. 



3. Strain-induced depolarization 
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The initial state is \m f i) and the final state is \m I j), 
where i and j represent nuclear spin eigenstates. m 
stands for the initial orbital state (QD s orbital). n qu 
is the Bose-Einstein distribution for phonons with mo- 
mentum q and phonon branch v at temperature T. 
We consider both phonon absorption and emission pro- 
cesses, depending on the direction of the total magnetic 
field. N q — fiq V for phonon emission processes, while 
N q = n qv + 1 for phonon absorption processes. fd e f{E z ) 
is obtained from the deformation potential term and 
fpiezo(E^) is obtained from the piezoelectric interaction. 
T ep is the transition amplitude for the phonon-assisted 
spin flip processes. Here the hyperfine interaction not 
only induces electron nuclear spin flip-flop, but also mixes 
electron spin and orbital degrees of freedom.!^ The elec- 
tron spin-up (-down) state in the m orbital is mixed 
with the electron spin-down (-up) state from a higher 
orbital state I. E m and Ei represent the energies of the 



Due to the strain-induced quadrupole interaction, 
where the principal axis Z is generally not parallel to the 
external field direction z, the nuclear spin eigenstates are 
a mixture of I z eigenstates. This means that the non flip- 
flop term in Eq. ©, I%S Z , can now induce transitions 
between different nuclear spin states and cause NSP (see 
red two-way arrows in Fig. [S]). Since the energy transfer 
between these nuclear spin states is generally much less 
than cotunneling energy (0.033/ieV^S), this process is not 
limited by energy conservation considerations. These de- 
polarization rates can be calculated by the Fermi Golden 
Rule, 



A As \ 2 2T c f e \(i\I z <-' 2 
N 



(7) 



where i and j represent nuclear spin eigenstates, and Ef 
and -E 1 ™ stand for the nuclear spin eigen- energies. 



D. Master equation of population 

Depending on the helicity of the excitation photon 
(<7 + /<7~), the electron can pump nuclear spins either to 
the higher-energy spin states or the lower-energy states. 
Take for example nuclear spins being pumped to the 
highest-energy nuclear spin state, as shown in Fig. [5l the 
average NSP can be evaluated by the master equation of 
population [see Eq. (0)]^ which is determined by the bal- 
ance between the pumping and depolarization channels. 



^Wi,jPi, 



(8) 
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where Wj t i ^ Wi.j represents the total transition rate where 
between the i and j states, and pi represents the nuclear 
spin population at ith state with i = 1, 2, 3, 4. 
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Here D'^ = Dij + is the total pumping rate. The 
system is highly nonlinear due to the population depen- 
dence of the pumping rate. The steady state nuclear 
spin polarization (I z ) — ^ li Pi (h) has to be calculated 
self-consistently. 
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IV. RESULTS OF NUCLEAR SPIN 
POLARIZATION 

Recent experiments have achieved NSP in In^Gai-^As 
QDs at various magnetic fields. NSP of 10 ~ 20% 
is created in Ino.9Gao.1As QDs at below 1 T^ 40% 
in Ino.6Gao.4As at approximately 2 X, 26 ' 27 and 80% 
in I1io.9Gao.1As at 5 T. 28 To better understand the 
differences in these results, we calculate the NSP 
in Ino.25Gao.75As (QD1), lno.6Gao.4As (QD2), and 
Ino.9Gao.1As (QD3) quantum dots at various magnetic 
fields. The default temperature for our calculations is 4 
K, unless otherwise identified. 

The z component strain tensor ezz is a good indica- 
tion of strain strength in our simplified model of the QD. 
As shown in Fig. [TJ the largest electric field gradient is 
Vzz, which is proportional to ezz- Near the surface of 
a pancake-like pure InAs QD embedded in GaAs, ezz 
ranges between 4 ~ 8% i 30 ' 36 i 37 ' 40 Inside, the strain is 
distributed more or less evenly, and ezz is 3.5%, which 
is about half of the lattice mismatch between InAs and 
GaAs (7%). To give a quantitative estimate, we assume 
that the strain information is completely contained in the 
ezz tensor element, which reflects the calculated lattice 
mismatch in dots with In^Gai^^As compositions^ 3 - By 
interpolation, we estimate ezz in QD1, QD2 and QD3 
to be 2.5%, 4.3%, and 6.3% respectively. 

It is worth noting that the electron g factor in a 
In-rGai-^As self-assembled quantum dot depends on 
the strain strength^ Based on the g factors given in 
experiments ) 26 ' 49 ^ 5 we take the g factors in QD1, QD2, 
and QD3 to be -0.6, -0.7, and -0.8 respectively. 



F 24 

Figure 6. (color online). Sketch of possible transitions be- 
tween nuclear spin states. The blue one-way arrows represent 
the pumping channels, and the red two-way arrows represent 
the depolarization channels. The thickness of the arrows qual- 
itatively represents the strength of the corresponding transi- 
tion. 

A. Single Arsenic nuclear spin polarization 

To understand nuclear spin polarization in different 
compositions of In^Gai-^As quantum dots, we first start 
with the nuclear spin polarization of As. The sim- 
plest case is the high magnetic field region, as shown 
in Fig. GDJd), where the nuclear spin Zeeman energy is 
much larger than the quadrupole splitting. Here nuclear 
spin eigenstates are close to the eigenstates of I z , so that 
{i\I z \j) w 0. Accordingly, the depolarization rates [see 
Eq. ([7])] between nuclear spin states are approximately 
zero. Therefore nuclear spins can be pumped to the high- 
est nuclear spin state, and nearly full nuclear spin polar- 
ization can be obtained. 

In low to intermediate field regions, the calculation for 
NSP becomes more complicated since the mixing between 
different nuclear spin Zeeman states become stronger 
than the higher-field case. The physical picture of various 
possible transitions is given in Fig. |51 

According to Fig. [SI in the absence of all the depolar- 
ization channels, nuclear spins can always be pumped to 
the highest spin state, and thus full NSP can be obtained. 
In the absence of -F12, no matter how strong other depo- 
larization channels (-F23, -F24, and -F34) are, the nuclear 
spins can still be pumped to state 1 eventually, and thus 
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Figure 7. (color online). An example of pumping and de- 
polarization rates in QD3. Fia and Fu do not appear here 
since they are smaller than 10 _5 /s. The state 2-4 degeneracy 
(where the nuclear Zeeman and quadrupole splittings coin- 
cide) occurs at approximately 0.67 T, where F24 and F23 reach 
their peaks. The pumping rates D + F (self-con) are calculated 
to self-consistency. In other words, these are pumping rates 
when nuclear polarization is already built up. D + r(random) 
are the initial pumping rates with a random distribution of 
nuclear spins before the Overhauser field builds up. 



become fully polarized. Once F12 is turned on, and in 
combination with F23 or F24, the pumped nuclear spins 
in state 1 can now leak back to state 3 or 4, and full 
polarization cannot be achieved. In other words, F12 is 
the key to depolarization. As an example, we plot all the 
depolarization rates and pumping rates in Fig. [7J 

In the two limits where E™ > Eqi and E™ < Eqi, 
the state spectrum is mostly Zeeman-like [Fig. [3Jd)] or 
quadrupole-like [Fig. G^a)]. In these cases F23 and F 2 4 
are very small compared to F12 (see Fig. [7]), so that pop- 
ulations pumped into states 1 and 2 cannot leak to states 
3 and 4. We can now simplify the four-level problem to 
a two-level problem, and find the steady state solution 
from the following equations: 



Pi 
Pi 



-F12 pi 
Vi = 1 ■ 



(D 1 2+T 1 2+F 12 ) P 2= 0, 



When E% < Eqi, the cotunneling- assisted spin flip tran- 
sition D\2 dominates, while when E" > Eqi the phonon- 
assisted spin flip transition is dominant, especially 
above 5 T (see Fig. Q|. At fields lower than approxi- 
mately 5 T, the average nuclear spin polarization can be 
expressed as follows, 



while at higher fields 



F12 
ri2+.Fi2 



(9) 



(10) 



From Eq. © and Eq. (J7J), when E™ > Eqi, F 12 is ap- 
proximately 
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and when E% < Eqi, 
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Equations ([9]) and (fT3|) can give a very good qualitative 
explanation to our calculations. At around 1.5 to 2 T, 
as shown in Fig. |8ja) and (d), D\% is at least one order 
of magnitude smaller than F\%, and the resulting (I z ) 
from Eq. © is - 1 (NSP is 67%), and is independent 
of field gradients of different angles and strengths (which 
determine -F12). 

There are some general trends in the NSP as evident 
in Eq. © (we focus on the regime of B < 5 T for our 
qualitative discussion in the following paragraph) , where 
(J 2 ) depends only on the ratio of depolarization to polar- 
ization F12/D12: 



1 



F12/-P12 

1 + F 12 /D 12 



For E» > Eqi, 
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and for E" < Eqi, 



F\2 

~D^2 



(f)2 + (-g*n B B z + 5 n f 



H\2 



(12) 



Notice that for _E™ > Eqi the ratio F12/D12 is propor- 
tional to the square of Eq/E™ and 8. When Eq or 9 
increases, -F12/-D12 increases, and the average NSP (I z ) 
will decrease. This is illustrated in the overall trends of 
Fig. [8j Likewise, in the regime of E™ < Eqi, the ra- 
tio F12/D12 is proportional to 9 2 and is a function of 
t c . Now when 9 increases, (I z ) decreases, again shown in 
Fig.[H] Furthermore, when the second term of Eq. (fl~2"|) is 
greater than the first (i.e. the electronic Zeeman energy 
is greater than the cotunneling energy), -F12/-D12 T c > 
so that (I z ) decreases when r c increases. As shown in 
Fig. GO the decrease of the cotunneling time constant (r c ) 
enhances the overlap of the electronic energy levels, and 
increases the cross-section of the hyperfine flip-flop pro- 
cesses. Therefore the resulting NSP increases, as shown 
in Fig. H 

The nonlinear nature of our system becomes most 
prominent when the Overhauser field is anti-parallel to 
the external magnetic field, as shown in Fig.[5Ja-c). Espe- 
cially when the Overhauser field cancels out the external 
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Figure 8. (color online). Nuclear spin polarization of As in three different In^Gai-^As quantum dots, QD1, QD2 and QD3. 
The strain strengths in the dots are approximately proportional to the In concentration x. The solid, dashed and dotted lines 
stand for the NSP in QDI, QD2 and QD3 with x — 0.25, 0.6 and 0.9 respectively. 9 represents the angle between the field 
gradient and the external magnetic field (the z-axis). Here 6 = 2°, 5° and 9° are considered. Panels (a), (b) and (c) show As 
nuclear spin polarizations at various strain conditions when the Overhauser field is anti-parallel to the external magnetic field. 
Panels (d), (e) and (f) show the parallel case. Overall, stronger strain or greater angle between the field gradient and growth 
direction suppresses the nuclear spin polarization. 



magnetic field (i.e. 2A(I z ) max » g*fi B B ext ) in D 12 , as 
shown for example in Fig. UOf b). At this point the spin 
pumping rate D\ 2 reaches its maximum, which leads to 
the peaks around 3 to 4 T in Fig. [5] The peak nuclear 
spin polarization is 



\ zlrnax ~ g 6 2 ( fo )2 ( JL )2 

[(^) 2 + (B?) 2 J 

~ 9 * llB ft 

~ 2 A v ext . 

Therefore, when Eq increases, the peak polarization 
(Iz) max decreases. The corresponding external field B ext 
decreases as well, but that relationship is more compli- 
cated because E™ also depends on B ext . 

Equation <j9j) generally fits well with our numerical 
calculations, except for the case when the nuclear Zee- 



man energy is equal to the quadrupole energy. When 
the nuclear Zeeman energy is nearly resonant with the 
quadrupole energy, the nuclear spin polarization is harder 
to build up because the degeneracy between states 2 and 
4 [see Fig. [3]Jb)] causes a maximum in the depolarization 
rate F 2 4, which is orders of magnitude larger than all 
other transition rates. The presence of this large tran- 
sition rate equalizes the populations of states 2 and 4: 
P2 = Pa- Furthermore, in combination with F12, it also 
prevents the nuclear spins from accumulating in the high- 
est nuclear spin state. By setting up the master equation 
of Eq. ((5]), we find the steady state solution from the 
following equations: 

F 12 pi = (D 12 + F 12 )p 2 , 

P3 »P4, 

Pi +P2 +P3 +P± = 1- 
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Figure 9. (color online). Field dependence of As nuclear 
spin polarization at different cotunneling rates. As shown 
in panel (a), the decrease of the cotunneling time constant 
(r c ), or the increase of the cotunneling energy, enhances the 
overlap of the electronic energy levels, specifically when the 
external magnetic field is below approximately 5 T. The cross- 
section of the hyperfine flip-flop processes is increased, and 
the resulting nuclear spin polarization increases accordingly. 
Panel (b) is a zoom-in of panel (a) near the 2-4 degeneracy. 



The average NSP is, 



^2 Pi 



where 
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The achievable NSP at the 2-4 degeneracy point depends 
on the relative direction of the Overhauser field and the 
external magnetic field. When they are parallel, -D12 is 
at least one order of magnitude smaller than F12, thus 
the NSP at this degeneracy is only a few percent, as 
shown in Fig. [SJd-f). When the fields are anti-parallel, 
D12 may become comparable to F12, and the resulting 
NSP strongly depends on the ratio of F12 to Di 2 , as 
shown in Fig.^a-c) and Fig. [TW a). If this ratio is closer 
to 1 (such as for QD1 and QD2), the resulting NSP could 
be as high as 20%, much higher than that is achievable in 
the parallel field case. On the other hand, if this ratio is 
far above 1 (in QD3) , the resulting NSP is suppressed, as 
shown in Fig. [Tola), while the field-dependence becomes 
abrupt. 

In order to gain more understanding of the highly non- 
linear behavior of DNSP in QD3, we examine the time 
evolution of D12 at three different energy detunings, 8\ , 
82 and ^3, away from the Zeeman-quadrupole resonance, 



Figure 10. (color online). The As nuclear spin polarization 
near 2-4 degeneracy [Panel (a)] and near the peak polariza- 
tion [Panel (b)[. The two panels are detailed views of Fig. [8j a) 
in these two regions. Near 2-4 degeneracy NSP depends on 
the ratio F12/-D12, as shown in Eq. ([13} . In QD1 and QD2, 
F12/D12 ~ 1, and the NSP is above 10% even at the de- 
generacy point. In QD3, F12/D12 3> 1, NSP is suppressed 
at the degeneracy and the field dependence of NSP becomes 
abrupt. The NSP peaks in panel (b) around 3 ~ 4 T occur 
when the Overhauser field cancels out the external mangetic 
field. These peaks shift to the lower magnetic field and lower 
nuclear polarization in stronger strain (see detail in text). 



as shown in Fig. [TTJ According to panel (b), when the 
pumping rate is high enough to overcome the depolariza- 
tion, the Overhauser field starts to build up. When the 
Overhauser field cancels out the external magnetic field, 
the pumping rate reaches a maximum in the time evolu- 
tion of the system [the spikes shown in Fig. [TTTb)]. The 
Overhauser field quickly exceeds the external magnetic 
field, and then this pumping rate falls off to a steady 
value, in a short time correlated to the hyperfine energy, 
cotunneling time constant and f e . For S3, the NSP build- 
up time is approximately 10 seconds, and for 62 it is 30 
seconds. For Si, the pumping never manages to over- 
come depolarization within our simulation time (2000 s), 
and the NSP is limited to a few percent. As the detun- 
ing S approaches zero, F24 gets closer to its maximum, 
which is orders of magnitude larger than all other rates. 
In this regime, it takes longer and longer time to build 
up the nuclear spin polarization, until it is practically 
impossible — beyond tens of seconds, nuclear spin relax- 
ation channels that we do not consider, such as dipolar 
induced spin diffusion and direct spin-lattice relaxation, 
would have to be included for a complete physical picture 
to emerge. 

When the nuclear Zeeman energy is equal to twice the 
quadrupole energy, nuclear spin states 3 and 4 are degen- 
erate, as shown in Fig. [3{c). However, the NSP around 
this degeneracy is not as strongly affected as those near 
the 2-4 degeneracy. In this field region, F 2 4 and F23 
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Figure 11. (color online). Time evolution of nuclear spin polarization rate [Panel (b)] at different energy detunings [defined in 
Panel (a)] in QD3 near the 2-4 degeneracy. The solid line at approximate 0.67 T in Panel (a) stands for the 2-4 degeneracy, 
where As nuclear Zeeman energy is resonant with the quadrupole energy. Si, 82 and S3 are 1.7 neV, 1.75 neV, and 1.8 neV from 
this degeneracy, respectively. Panel (b) shows that NSP build-up time for 83 is approximately 30 seconds and 10 seconds for 
82- For Si, within our simulation time of 2000 s, the nuclear spin polarization is not built up. Clearly, as 8 decreases, F24 gets 
closer to its maximum, which is orders of magnitude larger than all other rates, it takes longer and longer time to build up the 
NSP, until it cannot be built up. 



nearly vanish. The strong .F34 equalizes the population 
of states 3 and 4, so that we can again isolate the four 
level problem to a two level one, and NSP can be calcu- 
lated by using Eq. [§J 

B. Nuclear spin polarization in different 
compositions of In^Gai-a, As quantum dots 

In all the NSP experiments in In^Gai-^As QDs, the 
observable quantity is the total Overhauser field from all 
the nuclei of all the nuclear spin species. Since different 
nuclear isotopes have different gyromagnetic ratios (see 
Table [TJ) , and different nuclear spins generally experience 
different electric field gradients, the DNSP features we 
study in the previous section for a single nuclear spin 
species with a fixed quadrupole splitting would now oc- 
cur in ranges of magnetic fields. The total effect is a 
superposition of contributions from all the individual in- 
gredients. As we have discussed in section IIII1 we do 
not consider interactions between nuclear spins, whether 
they are of the same or different species. 

To account for the distribution of strain in a QD, we 
assume a uniform distribution of angles between the elec- 
tric field gradient and the applied magnetic field (which 
is along z direction as always). In Fig. the angles 
between the electric field gradient and the external mag- 
netic field are in the ranges of 0° - 5° and 0° - 9°. The 
peaks and dips in NSP of a single nuclear spin species, 
as shown in Fig. [5J are now smoothed out, as shown in 

pig.na 



Our results show a qualitative agreement with various 
experimentSi 26 i 27 ' 49 A high degree of nuclear spin polar- 
ization can be created in high field regions, while the 
polarization is limited in low field regions. Overall the 
achievable nuclear spin polarization in In^Gai-^As QDs 
is related to the concentration of indium and the result- 
ing strain distribution in the dots. In general, stronger 
strain and larger angle between the field gradient and 
growth direction suppress the nuclear spin polarization. 

V. DISCUSSION 

In our calculations of nuclear spin polarization, the 
spin pumping is achieved by first optically orienting the 
electron spins, then transferring the electron spin po- 
larization to the nuclear spins via either cotunneling- 
assitsed processes or phonon-assisted processes. In the 
latter we have not included any spin mixing caused by 
spin-orbit interaction. However, it is well known that 
the spin-orbit interaction is quite strong in InAs dotsi*^ 
Thus we have also explored how the spin-orbit interac- 
tion miff ht take part in the DNSP^l^ More specifically, 
we have calculated the spin-flip transition rate from the 
combination of the hyperfine interaction, spin-orbit in- 
teraction, and the electron-phonon interaction, and we 
find the transition rate is linearly proportional to the 
electronic Zeeman splitting. Our results show that the 
inclusion of spin-orbit interaction into the spin transfer 
process yields a slower process (by two orders of magni- 
tude even at the relatively low magnetic field of 1 T) than 
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Figure 12. (color online). The nuclear spin polarization in Ino.25Gao.75 As (solid lines), Ino.6Gao.4As (dashed lines) and 
Ino.9Gao.1As (dotted lines) quantum dots. Panels (a) and (c) are for dots with 6 6 — 5° while Panels (b) and (d) are for 
6 € — 9°. Panels (c) and (d) are zoom-ins to the low field region of Panels (a) and (b), respectively, (a) represents the 
nuclear spin polarization average over angles between 0° to 5° . (b) represents the nuclear spin polarization average over angles 
between0° to 9°. (c) and (d) are the detail in low fields for (a) and (b) respectively. Comparing scales of panels (a) with (b), or 
(c) with (d), we note that stronger strain and larger variations in the direction of the field gradient suppress the nuclear spin 
polarization in low to intermediate fields. 



the hyperfine interaction combined with electron-phonon 
interaction alone. Therefore, we exclude this mechanism 
from our current study. 

Our calculations presented so far are done at T = 4 
K. We have also explored the temperature dependence 
of the As NSP in a QD. Both cotunneling and phonon 
emission/absorption (especially absorption) are affected 
by temperature changes, via Fermi level broadening and 
phonon populations, so that spin-flip rate will change 
accordingly. Consider for example nuclear spins being 
pumped to the highest-energy state at high fields, where 
spin transitions are assisted by phonon absorption, as 
shown in Fig.[SJ We have calculated As NSP at three dif- 
ferent temperatures: 0.1 K, 4 K and 60 K. The results are 
shown in the high field region of Fig. Q21 where the NSP 
can be built up to larger values at temperature 60 K than 
4 K and 0.1 K (especially at B > 5 T), since at a fixed 
magnetic field the phonon population decreases as tem- 
perature decreases (N ~ cxp (g/fc B T)-i )' ^ n ^ e l° w field 
region, the spin transition is assisted by cotunneling, and 
the cotunneling time constant r c is inversely proportional 
to the temperature^ Therefore the cotunneling-assisted 
processes are more efficient at higher temperatures^ as 
shown in Fig. Q2J At present we do not have a clear an- 



alytical understanding of the abrupt change in NSP as 
shown in Fig. 1111 We are currently working on a full 
density matrix method that includes all the off-diagonal 
terms for the electron-nuclear spin system. Such a cal- 
culation could also help us ascertain the validity of the 
Master equation approach we adopt in the present study. 



VI. CONCLUSIONS 

In summary, we have calculated nuclear spin polar- 
ization through optical orientation of electron spins in 
a self-assembled quantum dot. We have explored how 
NSP of a single nuclear species depends on the exter- 
nal magnetic field with various strain strengths, angles 
between the electric field gradient and the growth direc- 
tion, and cotunneling energies. We show that, in high 
magnetic fields, higher degrees of NSP can be achieved, 
where the nuclear spin Zeeman energy is much larger 
than the quadrupole splitting. In this regime the nu- 
clear spin eigenstates are close to the eigenstates of I z , so 
that the depolarization rates between nuclear spin states 
are approximately zero. Therefore nuclear spins can be 
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Figure 13. (color online). Temperature dependence of As 
nuclear spin polarization. The solid, dashed and dotted lines 
stand for the nuclear spin polarization at temperatures of 0.1 
K, 4 K, and 60 K respectively. The nuclear spin polarization 
is easier to build up at higher temperature, since the phonon 
population and cotunneling constants increase as temperature 
increases. 

pumped to the highest nuclear energy state without leak- 
ing back to lower energy states. In low to intermediate 
field regions, NSP is strongly affected by the strain distri- 
bution. Generally speaking, in the same QD, the NSP is 
lower when the electric field gradient is at a larger angle 
from the external magnetic field, because strain along 
transverse directions (relative to the magnetic field) is 
the driving force behind depolarization transitions for the 
nuclear spins. In addition, NSP is lower at smaller cotun- 
nelling rates (when, for example, the allowed electronic 
state in the QD is far below the Fermi sea). Further- 



more, NSP is also harder to build up in a QD with a 
larger magnitude of strain. Our calculation shows that 
higher strain strength in a QD leads to smaller NSP in 
general. 

For NSP in In^Gai-^As quantum dots, our results 
are obtained from incoherent superpositions of In, Ga 
and As contributions in different proportions. We show 
that nearly full nuclear spin polarization can be created 
in high field regions, while it is limited in low field re- 
gions. Our results indicate that the concentration of in- 
dium and the resulting strain distribution in the dots 
play a crucial role in DNSP. For example, at low mag- 
netic fields, nuclear spin polarization is harder to build 
up in Ino.9Gao.1As than in Ino.6Gao.4As. The interplay 
between the nuclear quadrupole interaction and Zeeman 
splitting could lead to suppression of nuclear spin polar- 
ization. Our results are in qualitative agreement with the 
measured nuclear spin polarization in the experimental 
work of various groupsi 2 ^ 28 ! 32 ' 49 

Our results suggest that for a dot with a uniform strain 
distribution (and with a principal axis away from the 
external magnetic field), a minimum in NSP should be 
expected when the nuclear Zeeman energy is equal to the 
quadrupole energy. Moreover, a peak should be observed 
in the intermediate field regions (around 3 to 4 T), where 
the Overhauser field cancels out the external magnetic 
field. 
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